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Abstract 

Using the fractional derivative, considered in the Caputo sense, we study an 
analytical technique associated with the variational iteration method for the 
fractional generalized a-time Burgers’ equation with a > 0 and obtain approx¬ 
imate solutions in particular cases 0 < a < 1 and 1 < a < 2. 
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1. Introduction 

The fractional calculus (FC) is a very important tool associated with several 
problems which appear in physics, engineering an other sciences [I1I1I3I1]- 
As an important example we mention diffusion processes, particularly to the 
fractional partial Burgers’ equation (FPBE) appearing in the traffic flow and 
gas dynamics [S]. On the order hand, the variational iteration method (VIM) 
is a relatively new approaches to provide an analytical approximation to linear 
and nonlinear problems laiiiiH]. Those authors consider the VIM applied to 
the time-fractional partial Burgers’ equation. 
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Here we are interested in the VIM associated with the so-called generalized 
FPBE 


[w] — Uxx u — n(xjt), 


( 1 ) 


with D=^,0<a<l and l<a<2,HsM and p > 0. 

Some particular results appearing in the literature, are recovered. 

The paper is organized as follows: in section 2, preliminaries will be presented, 
particularly, a short review involving the FC, specifically the Caputo derivative 
and its respective Laplace transform, after of VIM it is analysed, in particular for 
nonlinear fractional partial differential equation. In section 3 we will present the 
VIM and the Burgers’ equation and provide a lemma, with the proof involving 
the general case and showing some particular cases of the parameters presented 
as examples. In the section 4 we will discuss the approximate solutions for 
the generalized FPBE, recovering the results involving the classical Burgers’ 
equation, presenting several examples with the respective graphics. Concluding 
remarks close the paper. 

2. Preliminaries 

In this section we will present definitions and results that we use in the 
paper, a short review of FC and the VIM for the linear and nonlinear equations. 

2.1. Fractional Calculus 

First of all, we introduce the Riemann-Liouville (RL) fractional integral, 
considered in the left, only j^]. Let n be a positive integer and a G C such that 
Re(a) > 0, we define the RL fractional integral by means of 



n — 1 < q; < n. (2) 
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The Caputo derivatives has been used by many authors in several physical 
applications [loiiniiniiiaiiiiiiiiiTe]. One reason for this choice is the fact 
that the initial conditions associated with the fractional differential equation 
are usually expressed in terms of integer order derivatives. Let n be a positive 
integer and a G C such that Re(a) > 0. We introduce the fractional derivative 
of order a in the Caputo sense, denoted by cD“[/(t)], by means of the integral 


1 


mim = 


_ r &Hr) 

r{n-a)Jo (t-T)“+i- 


-dr, n — 1 < a < n, 

I ' ' 


( 3 ) 


D"/(t) = a = n. 


The relation involving the Caputo derivatives and the RL fractional integral is 
given by 

,Dn/(i)] = ir“-cDr/(t) (4) 


with the n — 1 < a < n. 

As we have already said, the Laplace transform methodology is an efficient tool 
to discuss a fractional differential equation. Then, we introduce the Laplace 
transform of the derivative in the Caputo sense. Denoting by ^ the Laplace 
integral operator, we can write the Laplace transform of the Caputo derivatives 
as follows 

n—1 

(.5-[,D“/(t)])(s) = s“(if[/(t)])(s) - ^ s“-'=-1(DV)(0+) , (5) 

fc=0 

with s is the parameter of the Laplace transform. 

As we have mentioned above, this expression shows that the Laplace transform 
of the fractional derivative in the Caputo sense involves only the derivative of in¬ 
teger order evaluated in t = O'*', conversely the corresponding RL derivative. In 
our particular problems, as we will be seen in the sequence, we take the param¬ 
eter a as a real number such that 0 < a < 1 in problems involving (anomalous) 
diffusion and 1 < a < 2 in problems associated with wave propagation. 

To close this subsection, as an example, we consider the fractional integral and 
the fractional derivative of a power function with A a real parameter. For 
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the fractional integral we have 


.ar.Ai ^ r(A + 1) 

^ r(A + l + a) 


( 6 ) 


with a > 0, A > —1 and t > 0. On the other hand, for the fractional derivative 
we have 




r(A + i) 
r(A + l-a) 


(7) 


with a > 0, A > —1 and t > 0. 


2.2. Variational iteration method 

The VIM [TTKISlIin] was extended to fractional differential equations and has 
been one of the methods frequently used. Classical and fractional differential 
equations are studied using VIM [181 HI]- On the order hand classical and 
fractional partial differential equations are studied in ISQlEllEaEaEllEgEel 
HZIIIH] in particular nonlinear dynamics for local fractional Burgers’ equation 
arising in fractal flow is discussed in [2^ . 

Here we consider a more general fractional differential equation 

cDf[u] + i?[M] + N[u] = f{t), 

where cD“[m] is the Caputo derivative, is a linear term, N[u] is a nonlinear 
one and f{t) is a function associated with the non homogeneous term. Odibat 
and Momami in [301 applied the VIM to the above equation and suggested a 
variational iteration formula 


Un+i =Un+ [ A(t, r)(cD“M„ + R[u] + V[m] - /(T))dr, 0 < t 

Jo 

' X{t, t) = —1, 0 < a < 1 

A(t, t) = T — t, 1 < a < 2. 

A(t, t) are known as the Lagrange multipliers associated with the variational 
iteration formula, this multipliers are evaluated with the general theory of La¬ 
grange multipliers [31j . 
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3. VIM and fractional Burgers’ equation 


In this section, we will present the VIM applied to the Burgers’ equation: 


,D“[u] = Uxx + AvFux, u = u(x, t) 


(8) 


with D = —,0<a<l and l<a<2, AgK and p > 0. 
at 


The so-called correction functional for Eq.(8) is 


Uk+l{x,t) = Uk{x,t) + 


r(/3) 


{t-rf 1 \{t) ( ^Uk - ^Uk - Aul^uk ) dr, 


A 

‘ dx' 


(9) 


where j3 = a + 1 — m and m — \ < a < m. 

The Eq.(9) can be approximately expressed by means of 

ft / Qra q2 ^ _ d _\ 

Uk+i{x,t) = Uk{x,t) + j X{t)\—-^Uk--f^Uk-Au\—Uk\dT. (10) 

If we use integration by parts and remembering that the stationary term in the 
functional 5uk = 0, we get three cases as follow: 


a) For TO = 1 we have 0 < a < 1 and the correction functional can be approxi¬ 
mately by means of 

duk+iix,t) = 6ukix,t) + S J A(t) ( — Ufc - Mfc 1 dr. 

Thus, we have the system 


whose solution is A(r) 


1 -H A(r) = 0 
A'(r) = 0, 


( 11 ) 


— 1. Then, we obtain the following formula (note that. 
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here we have j3 = a) 


If* / ,9“ 92 d \ 

Uk+i{x,t)=Uk{x,t)- ^ g^**k-Aul—Ukj dr. 

( 12 ) 

b) For m = 2 we have 1 < a < 2 and the correction functional can be approxi¬ 
mately expressed by means of 


6uk+i{x,t) = 6uk{x,t) + 6 X{t) 


52 g 


dr. 


Thus, we get the system: 


A"(r) = 0, 

XiT)r=t = 0, (13) 

1 - A'(r),=i = 0, 

whose solution can be written as A(r) = t — t. Then, for m = 2 and 1 < a < 2, 
we obtain the following interaction formula with (3 = a — 1 


, , , , a — lf*. / 5 “ A V ^ \ ^ 

Uk+i{x,t) = Uk{x,t) - [gt^^k --^Uk - Au'f—Uk j dr. 

(14) 

c) News Lagrange multipliers with 0 < t and 0 < a. In this case, our main 
result, we propose a lemma with its proof, recover the two precedent results 
and present an example. 


Lema 1. If the correction functional of the Eq.(9) is given by the Riemann 
integration 


Un+l=Un+ X{t,T)[cDf{Un)-{Un)xx-A{Uny{Un)x]dT, (15) 
Jo 
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with 0 < t, 0 < a and (un)xx, A(unY(un)x are restrictions of the variations of 
the functional associated with Eq.(15), then the Lagrange multiplier is 


X{t,T) 




Proof. First of all, we transform Eq.(15) in its integral form and taking the 
Laplace transform on both sides of the new equation 


^[Un+i] = ^[Un + \ fXRL{t, t)(cD“u„ - {Un)xx “ A{UnY{Un)x)] (16) 


where XrrY^t) is the Lagrange multiplier of the integral form for Eq.(15). We 
consider 


\fXRL{t,T)cDfUn= [ {t-T)° ^A(t, T)cD>„(T)dT. (17) 

r(a) Jo 


Thus, if XrlYjT) = X{t — r), Eq.(17) is the convolution of the function 

A(t)t“-i 


i{t) = 


r(a) 


(18) 


and cD“u„(t). The terms {un)xx and A{unYiun)x which are considered restric¬ 
tions on variations, implying S{un)xx = 0 and dA{unY{un)x = 0. 

Using the variational functional associated with Eq.(16) and the Laplace trans¬ 
form of the Caputo derivative, we obtain 


SJ^[Un+l] = S^[uY + [l“AflL(cD“u„ 


Yn)xx Al(rt^)^('Uyj)x)], 


m —1 

= 6J^[Un] + 5[J^[a{s)]s°‘J^[Un{s)]\ — 5 ^ u*(0''')s““^“^ 

k=0 


= (l+^[a(s)]s“)<5^M, 


(19) 
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where, using the Euler-Lagrange equation associated with Eq.(16), we have 


1 + ^[a(s)]s“ = 0, 

resulting 

( 20 ) 

Performing the inverse Laplace transform of Eq.(20) we get 

= ( 21 ) 

r(a) 

Comparing Eq.(18) and Eq.(21) we have r) = —1 and using Eq.(16) we 

can write 


'^n+1 — '^n 


Un-\-l — '^n Ufi (Un^xx ('Wn)a:] ? 

/•* (t-T)“-l 


^n +1 — 


Jo r(a) 

(-l)“(r-t)“-i 

TH 


Un {Un)xx Al(ri^)^(Wyj)a,]dT, 


(Wyi) {Un)xx Al('Un)^('U7^)a;] dT. 


Finally, using Eq.(17) we obtain the following general expression 

(-l)“(r-t)“-i 


A(t,r) = 


r(a) 


( 22 ) 


In what follows, we recover the known results and discuss an example. 


Corolario 1. In Eq.(22) if a = 1 then A(t,r) = —1 and if a = 2 then 
X{t, t) = T — t. 


Proof. Direct substitution a = 1 (a = 2) in Eq.(22). 



Example 1 : Analogously to the results that can be obtained by [8], we have the 
approximated solution for the FPBE when 0 < a < 1 using the initial condition 
u{x, 0) = g{x). The problem to be consider is 


,Dfu + u 


du 

dx 


d^u 

= 0, t > 0, 0 < a < 1, 


(23) 


[ u{x, 0 ) = g{x), 0 < a; < 1. 

Using Corollary 1 we have X{t,T) = — 1 and substitution in the correction 
functional Eq.(15) we obtain 


Un+l{x,t) 




O-trv ^n{Un)a 

ot°‘ 



The correction functional can be written as 

rt 


1 f 

'^n+l (^7 0 ~ 0 pO C / T '^n ('^n)x ('^n)a:x]dT. (24) 

r(a) Jo 

Using Eq.(24) we obtain the following approximations 


uo = u(x,0) = g(x), 


ui = g{x) - [gg' - 5 "] 


r(a + l)’ 


U 2 =u, + {2g{g'f + g^g" - 2gg^^^ - 4g'g" + 


t 


.2a 


r(l + 2a) 


-t3a 


Similar calculus are using for obtain others approximations. To close the section 
we prove another corollary. 


Corolario 2. With \{t,T) - 
Un+l {x,t) = Un {x,t) + 


(-l)“(r-t)“-i 

TM 


and correction functional as 


I (t-r)^ ^X{t,T)[cD“-{Un)xx-A{UnyiUn)x]dT 
0 

(25) 
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where (3 = a — [cx\ and [a] is integer part of a, then the iteration equation for 
0 < a < 1 and 1 < a < 2, is 


'^n+l — '^n If [cl^f id^n)xx 5 

and 

'^n+l — (Q^ ^)lf {c^t {Un'jxx (^Un) x] ■ 

respectively. 

Proof. By the Corollary 1, \{t, r) = —1 for 0 < a < 1 therefore jd = a, which 
imply 


1 /■* 

Un+l{x,t) =Un{x,t) - / (t - r)“"^ [cD“m„ - yl('U„)^’('U„)^ 

r(a) Jo 

or in the following form 


i^n) xx\^'^ 1 
(28) 


'^n+l{x.ft') — Uji 1^ [cD^Uyi {Un)xx ('^n)x] ■ 

Now for 1 < a < 2 therefore /3 = a — 1, by the Corollary 1, we have \{t, r) = 
T — t and using the relation 


1 -(a-1) 

r(a- 1) “ r(a) ’ 

we obtain 

1 /■* 

Un+l {x^t') = Uji / {t t) [cDf {d^n)xx {’^n)x]^X , 

r(a- 1 ) Jq 

or in the following form 

^n+l {x^t) = Uji (o 1)1^ [cDf id^njxx id^n)x\ • (^0) 
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4. Approximate solutions for FBPE 


There are some researchers that consider fractional Burgers’ equation to model 
the diffusion behaviour of the flow through porous medium. In this section we 
will consider three examples of the fractional Burgers’ equation, two of them 
with 0 < a < 1 and another one for the case 1 < a < 2 , with its respective ini¬ 
tial conditions. Here u is the flow of velocity, the viscosity coefficient is consider 
equal to —1, and without loss of generality we take A = —1 and p = 1 in the 
Eq.(l). Note that, the viscosity coefficient corresponds the second derivative in 
the Eq.(l). 

For those three examples, we consider an approximation of the solution in the 
following series form 

u{x, t) = uo{x, t) + ui{x,t) + U 2 {x,t) + ■ • • 

To make the graphics, we stop the series in the three term. 

We mention that the effect of the fractional derivative recover the memory 
effect associated with physical phenomena for 0 <a<l,l<a <2 [5?| . 
In particular with suitable initial conditions, when A = 0, a = 1 and a = 2 
in Eq.(l), we recover the memory effect associated with a heat equation and 
wave equation, respectively. Also, with a = 1 in Eq.(l), we obtain the classical 
Burgers’ equation. 

To close this section we will present some examples with graphics to see better 
the effect involving the parameter associated with the derivative. 


Example 2: Consider the problem 


r.™ du d'^u 


l{x, 0) = g{x) = sin(7ra;). 


t > 0, 0 < a < 1, 

0 < a; < 1. 


(31) 
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Using Eq.(25) with A(t,r) = —1 we obtain 


uo(x, t) = sin(7ra;), 


ui(x,t) 


sin(7rx) — 7rsin(7rx)[cos(7ra;) + (tt)^ 


r(a + l)^ 


U 2 (x,t) 


ui(x,t) + (27r^ sin(7ra;) cos^ (ttx) 


TT^ sin^(7ra;) + 27r^ sin(7ra;) cos(7ra;) 


+47r^ cos(7ra;) sin(7ra:) + 


4 ■ 4 

TT sin 


(ttx)) 


r(l + 2a) 


—7rsin(7ra;) cos(7ra:) + (tt)^ cos(7ra;) — tt^ sin^(7ra:) 


—TT^ sin(7rx)) 


r(l + 2a)t3“ \ 

r2(i + a)r(i + 3a)y' 


The graphics in Figure 1, with a = 0.2 and Figure 2 with a = 1 elucidate the 
approximation of u(x, t). 
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Figure 1: u(ir, 0) = sin(7rx), with a = 0.2. 



Figure 2: ii(ic,0) = sin(7ra:), with a = 1. 


Example 3: Consider the Example 2 with the initial condition given by u{x, 0) = 

g(x) = -. Using Eq.(25) with A(t,r) = —1 we get 

1 + 


uoix,t) = 


1 + e“’ 
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Ui{x,t) 


1 + 


V (l + e-)3 ) lvr(a + l)j’ 


U2{x,t) 


l + e“ 


/e"^(e“+ e2^ - 1)\ / ^ \ 

V (l + e“)3 j lvr(a + l); 


/e"=(l - 12e“ - 29e2^ + 4e3“ + 37e*^ + 50e^^ + lle^^ + 2e^^) 

(l + e^)8 ’ 


Vr(l + 2a) 


e^(l _ 4e^ _ I0e2^ + 

(1 + e®)3 


/r(l + 2a) t 3 a \ 

^ + a) T{1 + 3a)) ' 

Also here, the graphics in Figure 3, with a = 0.2 and Figure 4, with a = 1, 
elucidate the approximation of u{x,t). 



Figure 3: u[x^ 0) 


-, with Q = 0.2. 

1 + e^ 
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Figure 4: u(x, 0) 


-, with 0 = 1. 

1 + 


Example 4: Consider the problem 


du d‘^'11 

+ w—- — = 0 , t > 0 , 1 < a < 2 , 


M(a:, 0) = sin(7ra;), Mt(a;,0)=0, 0 < a; < 1. 


In this case, using Eq.(25) with A(t, t) = t — t we obtain 


uo(x, t) = sin(7ra:), 

ui{x,t) = sin(7ra;) — a{a — l)[7rsin(7ra;) cos(7ra;) — sin(a:)] 
U 2 (a:, t) = sin(7ra:) — a(a — l)[7rsin(7ra;) cos(7rx) — tt^ sin(x)] 




,r(a + i). 
r 

.r(a + i). 

+a{a— l)[sin(7rx)(27r^ cos^(7rx) — tt^ sin^( ttx) + Gtt^ cos(7rx) 
+7r^ sin^(7rx)] 


_r(l + 2a) 

X (tt^ cos(7rx) — TT^ sin^(7rx) — sin(7rx) 


— a{a — l)(7rsin(7rx))(cos(7rx) + x^) 

r(l + 2a)t^a 


_r 2 (i + a)r(i + 3a)_ ■ 
The graphics in Figure 5, with a = 1.2 and Figure 6, with a = 2 one can see 
the evolution of the approximation of u{x,t). 
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Figure 5: u(x, 0) = sin(7r3;), Mt(a;, 0) = 0, with a = 1.2. 



Figure 6: u{x, 0) = sin(7rx), ut{x, 0) = 0, with a = 2. 


5. Conclusions 

The FC is very useful in the recuperation of the memory of phenomena, by 
used of fractional derivative in time variable. News Lagrange multipliers for 
FPBE are identified by Laplace transform and particular cases are recovered. 
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Using this multipliers and the VIM we obtained approximations of the solutions 
for FPBE taking three terms only. Then, we conclude that the VIM is a powerful 
and efhcient technical to approximate solutions to FPBE. 
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